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High resolution heat capacity and thermal expansion experiments performed on single crystalline 
kagome-staircase compounds (Coi_a:Nii,)3V208 are presented. The parent compounds C03V2O8 
and Ni3V208 undergo a complex sequence of first- and second-order magnetic phase transitions. 
The low-temperature (T < 40 K) magnetic entropy evolves monotonously with the doping content x, 
from the full S = 1 Ni^""" magnetic entropy in Ni3V208 to half of the S = 3/2 Co'^"'' magnetic entropy 
in C03V2O8. Thermal expansion coefficients ai {i = a, b and c) show a strong anisotropy for all 
(Coi_a;Nii,)3V208 compounds. The low-temperature magnetic distortion indicates that Co-doping 
(Ni-doping) has similar effects to applying a uniaxial pressures along a or fe (c). Linear Griineisen 
parameters Ti are extracted for the three main axes i and exhibit a complex behavior with both 
temperature and doping. For each axis, Ti and a^ exhibit a sign change (at low temperature) at 
the critical concentration Xc — 0.25, at which the incommensurate magnetic propagation vector 
changes. Beyond our study, an understanding of the multiple and complex parameters (magnetic 
frustration, magnetic anisotropy, mixture of 5 = 1 and 5* — 3/2 ions, etc.) is now necessarily to 
bring light to the rich magneto-elastic properties of (Coi_i:Nia;)3V208. 
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PACS numbers: 75.30.-m, 75.30.Kz, 75.50.Ee, 65.40.-b 



INTRODUCTION 



The study of the kagome-staircase compounds 
C03V2O8 and Ni3V2 08, which have the same orthorhom- 
bic (Cmce) crystalline structure^i has recently stimulated 
both experimental and theoretical works. In these sys- 
tems, two-dimensional and geometrically-frustrated mag- 
netism results from a kagome-staircase arrangement (per- 
pendicularly to the 6-axis) of the magnetic ions. The 
magnetic lattice is composed of 5 = 3/2 Co^+ ions 
in C03V2O8 and 5 = 1 Ni2+ ions in Ni3V208. It is 
made of two non-equivalent magnetic sites. ^ In both com- 
pounds, cascades of magnetic transitions have been re- 
ported at low temperatureSj^ii. which indicates a high 
number of competing magnetic interactions. Both sys- 
tems are characterized by a "high temperature" in- 
commensurate (HTI) magnetic order (below Thti.i = 
11.3 K in C03V2O8 and below Thti,2 = 9.1 K in 
Ni3V208 [J]) and by a low temperature commensurate 
phase (below TpM = 6.2 K in C03V2O8 and be- 
low Tc = 2.2 K in Ni3V208 d). They also have 
the same easy magnetic axis a. However, their mag- 
netic ground states are associated with different order- 
ing wavevectors. For example, C03V2O8 is a weak fer- 
romagnet below Tpi^i — 6.2 K |3|, while Ni3V2 08 is 
a canted antiferromagnet below Tc' = 2.2 K [J]. In- 
terest for these systems was raised with the discovery 
that Ni3V208 enters into a multiferroic regime for Tc — 



3.9 K < T < Tlti,2 = 6.3 K, where canted antiferromag- 
netism is accompanied by a ferroelectric polarizationi^ 
Contrary to Ni3V208, C03V2O8 is not multiferroic 1^ Re- 
cent neutron^''* and specific heat^ studies permitted to 
draw the doping-temperature magnetic phase diagram 
of the mixed compounds (Coi_2:Ni^)3V208 and to de- 
termine the evolution of the magnetic structure of the 
system, from C03V2O8- to Ni3V2 08-type, as a function 
of doping. From these measurements, only one transi- 
tion was reported for the samples of nickel concentrations 
0.14 < a; < 0.73. A change of the magnetic structure was 
found to occur around the concentration Xc — 0.25 at 
which the transition temperature Thti — 5.5 K is the 
smallest,^ 

Here, we present a study by specific heat and ther- 
mal expansion of single crystals of the mixed compounds 
(Coi_j;Nij:)3V208 with nickel concentrations a; = 0, 0.12, 
0.2, 0.28, 0.35, 0.5, and 1. The same sequences of tran- 
sitions as those reported in the literature^Ti^i^iiSrJ^ are 
found for the parent compounds C03V2O8 and Ni3V208. 
The high resolution of our specific heat and thermal ex- 
pansion experiments permits us to study the temperature 
and doping evolution of these transitions in detail and 
to derive a detailed temperature-doping phase diagram. 
Further, magnetic Griineisen parameters, which provide 
a sensitive probe of the volume/pressure dependence of 
the magnetic energy scales, can be derived from the com- 
bination of thermal expansion and heat capacity data. 



We will demonstrate that the behavior of the Griineisen 
parameters in this system is extremely complex, proba- 
bly due to a high number of competing magnetic energy 
scales. It exhibits some parallels to the novel behavior 
expected at quantum phase transitions/'^ which will be 
discussed in detail. This paper is organized as follows. 
In Section [2] we provide experimental details. In Section 
|3] we present the specific heat and thermal expansion re- 
sults, as well as the calculated uniaxial Griineisen pa- 
rameters. In Section m the phase diagram is derived from 
our data and the pressure dependences and Griineisen 
parameters are discussed, and finally the conclusions are 
presented in Section [5] 



2. EXPERIMENTAL DETAILS 

The (Coi_a;Ni2;)3V208 single crystals with Ni-contents 
X = 0, 0.12, 0.2, 0.28, 0.35, and 0.5 studied here were 
grown from fiuxes composed of C03O4, NiO and V2O5. 
Ni3V208 single crystals were grown from a flux composed 
of BaO, NiO and V2O5. Ni-contents x were determined 
assuming the validity of the Vegards law, i.e., that the 
lattice constants (measured using powder x-ray diffrac- 
tion) vary linearly with x. The crystals have been cut 
to typical dimensions of about 5x5x2 mm'^. Specific 
heat was measured using a Physical Properties Measure- 
ment System (PPMS) from Quantum Design; a "^He in- 
sert was used to reach the lowest temperature of 0.3 K. 
Since the standard PPMS software is not adapted to the 
study of sharp transitions, a relaxation method similar 
to the ones proposed in Ref. [ 14" 15*1 was used to obtain 
reliable data near the first- and second-order phase tran- 
sitions. Linear thermal expansion was measured between 
2 K and 300 K usinga home-made high-resolution capac- 
itance dilatometer [16lll7| with heating and cooling rates 
of 6 mK/s or 20 mK/s, within a resolution AL/L = 10~^. 
For all samples, thermal expansion was measured along 
the three main crystallographic directions a, b, and c. 



3. RESULTS 

3.1. Specific heat 

Specific heat Cp versus temperature T data measured 
on C03V2O8, (Coo.5Nio.5)3V208, and Ni3V208 are pre- 
sented in Fig. [U for 0.3 < T < 30 K (log-log scale). 
The same sequence of transitions as those reported in 
the literature^Ti^'^ are observed for the two parent com- 
pounds C03V2O8 and Ni3V208. Huge first-order anoma- 
lies are obtained at Tp]\f — 6.00 ± 0.05 K in C03V2O8 
and Tc = 3.90 ± 0.05 K in Ni3V2 08, at which Cp 
reaches 285 Jmol^^K^^ and 64 Jmol~^K~^, respectively 
(almost an order of magnitude higher than in Ref. [2- 
|3la|)- The high quality of our single crystals and an 
appropriate procedure to analyze the data (relaxation 
methodiiiii) permitted to get these high-resolution data. 



(Coo.5Nio.5)3V208 has only one second-order phase tran- 
sition at 7.45 ± 0.50 K (see Appendix [X| . The specific 
heat follows a T^ power law below 1.5 K in Ni3V208 and 
a T'* power law for 2 < T < 4 K in C03V2O8, leveling off 
below 2 K. These laws presumably result from low-energy 
magnetic excitations. The T^ law in Ni3V208 is compat- 
ible with three-dimensional ungapped antiferromagnetic 
spin waves, in relation with the canted antiferromagnetic 
state which sets in below Tq ~ 2 K. Oppositely, the 
T^ behavior in C03V2O8 might result from more com- 
plex spin dynamicsi^ The cascade of magnetic transi- 
tions in C03V2O8 and Ni3V208 are also accompanied by 
an additional signal proportional to T. This linear Cp 
"background" might be related to additional spin fluc- 
tuations, similarly to what is observed in heavy-fermion 
antiferromagnets (see for an example Ref. [19u20] ). For 
the doped (Coo.5Nio.5)3V208 compound, the low tem- 
perature specific heat is enhanced, indicating stronger 
low-energy magnetic excitations than in the parent com- 
pounds. 

Fig. m shows a plot of Cp up to 300 K for C03V2O8, 
(Coo.5Nio.5)3V208, and Ni3V208, together with esti- 
mates of the phonon background. The non-magnetic 
signal was modeled by a combination of one Debye and 
two Einstein phonon modes, which were fit to the data 
above 60 K (within the approximation that the sig- 
nal is purely phononic above 60 K). The fitted phonon 
backgrounds were subtracted from the original data to 
yield the magnetic heat capacity Cmag shown in Fig. 
13] (a) as Cmag/T versus T, for 0.3 < T < 60 K 
((Coi_2;Ni^)3V208 samples of concentrations x — 0, 
0.12, 0.2, 0.28, 0.35, 0.5, and 1). From our specific 
heat data, C03V2O8 undergoes the series of transitions 
Thti,i = 11-26 ± 0.02 K (" high- temperature incom- 
mensurate" phase), Thtc,i = 8.50 ± 0.05 K ("high- 
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FIG. 1: (Color online). Specific heat of C03V2O8, 
(Coo.5Nio.5)3V208 and Ni3V208 as a function of temperature 
in a log-log scale. 



temperature commensurate" phase), Tlti,2 = 6.20 ± 
0.05 K ("low-temperature incommensurate" phase), and 
TpM — 6.00 ± 0.05 K (low-temperature ferromagnetic 
phase), while Ni3V2 08 is characterized by the transitions 
Thti,2 = 9.42±0.02 K ("high-temperature incommensu- 
rate" phase), Tlti,2 = 6.48 ± 0.02 K ("low-temperature 
incommensurate" phase), Tc,2 = 3.93 ± 0.02 K (canted 
antiferromagnetic phase), and Tc,2 — 2.42 ± 0.03 K 
(second canted antiferromagnetic phase) - see Appendix 
\M Thtcs, Tlti,2, and TpM in C03V2O8 and Tc,2 in 
Ni3V20g are of first-order 1^ In contrast to the end com- 
pounds, the heat capacities of the mixed crystals with 
X = 0.12, 0.2, 0.28, 0.35, and 0.5 show only one phase 
transition, which is of second order, at 8.19 ± 0.10 K, 
7.6 ± 0.2 K, 6.4 ± 0.4, 7.0 ± 0.4 K, and 7.9 ± 0.3 K, re- 
spectively. The transition temperature goes through a 
minimum at x — 0.28 %. For the mixed compounds, 
(jmag ^rp jncj-eascs with decreasing temperature below 
0.5 K. Stronger low-temperature magnetic fiuctuations 
might result from a stronger magnetic frustration and 
could explain the high value of the specific heat at 
low temperature. Similarly, enhancement of the low- 
temperature specific heat is generally observed at the 
quantum phase transition of heavy-fermion systems ^1^ 
where low-energy fluctuations are known to play a central 
role.22 

The magnetic entropy Smag (from an integration of 
(jmag jrp ygj.g^g j^) jg g^own for different x in Fig. [3] 
(b) . Clear step-like increases of the magnetic entropy are 
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FIG. 2: (Color online). Specific heat - and its non- 
magnetic contribution - versus temperature of C03V2O8, 
(Coo.5Nio.5)3V208 and Ni3V208. 
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FIG. 3: (Color online), (a) Magnetic heat capacity divided by 
temperature C^"'^ /T and (b) magnetic entropy versus tem- 
perature of (Coi_^Ni2,)3V208 single crystals. 



Co vo. 

3 2 8 



(Co, Ni ) V O 

^ 1-x x'3 2 8 



35 



^-~ 30 



15 



3Rln4 = 34.55 Jmol K (full entropy from S = 3/2 Co Ions) 



3Rln3 = 27.38 Jmol'^K"' (full entropy from S = 1 Ni^* Ions) 




0.0 



0.2 



0.4 0.6 

Ni content x 



0.8 



1.0 



FIG. 4: (Color online). Plot of the magnetic entropy of 
(Coi-2,Nia;)3V208, at T = 40 K, as a function of the Ni con- 
tent X. 



observed at the first-order transitions TpM in C03V2O8 
and Tc.2 in Ni3V208. In Fig. 21 the saturated mag- 
netic entropy (at 40 K, from Fig. [3] (b)) is plotted as a 
function of x. For Ni3V2 08, Smag equals at 40 K the 
full entropy 3i?ln3 from S = \ Ni^"*" ions, as already 
reported in Ref. @|j|- For C03V2O8, Smag reaches at 
40 K only half of the full entropy 3i?ln4 expected for the 
S = 3/2 Co^+ ions, which is also in accordance with pre- 



vious studies/^— Rogado et al.^ initially proposed that 
the missing entropy could be related to additional mag- 
netic transitions occurring below 4 K (4 K was the low- 
est temperature investigated in their study). We have 
shown here that the specific heat of C03V2O8 decreases 
strongly with decreasing temperature (down to 0.3 K) 
and the shape of our specific heat data indicates that 
no further transition might occur below 0.3 K. Alterna- 
tively, Yasui et al.^ proposed that a crystal-field split- 
ting of the S = 3/2 quadruplet leads to a ground state 
composed of an effective 5*6// = 1/2 doublet, associated 
thus with a magnetic entropy of 3-Rln2 (half of 3i?ln4). 
The hypothesis of a stronger crystal field in C03V2O8 
than in Ni3V208 is compatible with former susceptibil- 
ity measurements,—"— which showed that the magnetic 
properties of C03V2O8 are very anisotropic while those 
of Ni3V208 are almost isotropic. We note that a similar 
situation occurs in the chlorides C0CI2 and NiCl2 , where 
crystal-field effects are stronger in C0CI2 than in NiCl2 



3.2. Thermal expansion and Griineisen parameter 

Fig. [5] (a) and (b) show the temperature dependence 
of the thermal expansion coefficients a^ of C03V2O8 
and Ni3V208, respectively, for i — a, b, and c and 
2 < T < 300 K. They are defined as a^ = l/L,dL,/dT, 
where L,; is the length along i. The same sequences 
of phase transitions as those reported by specific heat 
(Section [2]) can be seen in our thermal expansion data 
for both C03V2O8 and Ni3V2 08 (see Fig. [3 Table U 
and Appendix |X]) . The thermal expansion anomalies in 
C03V2O8 are positive and of roughly the same magnitude 
for the a- and 6-axes and are negative for the c-axis. In- 
terestingly, this situation is reversed for Ni3V2 08, which 
has mostly negative effects along the a- and &-axes and 
positive anomalies along the c-axis. Above about 50 K, 
the magnetic contributions vanish and both compounds 
exhibit very similar expansion coefficients due to the an- 
harmonicity of the phonons. Here, ab and ac are nearly 
identical for both compounds and aa is somewhat larger. 
The similarity of the phonon thermal expansion for both 
compounds is not unexpected since Co and Ni are neigh- 
bors in the periodic table. 

The detailed temperature variation of the thermal 
expansion coefficients, the integrated relative length 
changes, as well as uniaxial Griineisen parameters are 
shown in Fig. [5] for various Ni concentrations x. We 
follow Ref. J^ in defining the uniaxial Griineisen param- 
eters, Fi, as just the ratio Ui/Cp. The normalized lengths 
[L, - L,(30 K)/L,(30 K) for C03V2O8 and Ni3V208 
and i = a,b,c again nicely demonstrate the anisotropy 
discussed earlier and are in good agreement with pre- 
viously reported values . ^^'^^ For the intermediate com- 
pounds, there is a progressive change from the Co-type 
anisotropic behavior to the Ni-type behavior as x in- 
creases. The mixed compounds all exhibit only the high- 
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FIG. 5: (Color online). Thermal expansion coefficients aa, ctb, 
and Qc versus temperature of (a) C03V2O8 and (b) Ni3V208. 



est temperature phase transition, except for the x = 0.12 
Ni-sample, which shows two transitions. The anomalies 
in tti become much smaller for the intermediate com- 
pounds, and the change of behavior from Ni- to Co-like 
occurs around x = 0.20-0.28. At these concentrations, 
the exact behavior is in fact quite strange; above Thti 
the data appear to belong to the Co-side, then at Thti 
and below Thti the a^ jumps and sign changes to the Ni- 
like behavior. This is most apparent for the a- and c-axis 
data (Fig. [6]d and f). We note that this 'jump' occurs 
only in the thermal expansion data and not in the heat 
capacity data, where the transition is only marked by a 
kink in Cp, and that, therefore, it is not at all apparent 
how one can apply the Ehrenfest relation for calculating 
the pressure dependences of Thti ■ 

In order to gain a better understanding of the above 
behavior, it is informative to examine the behavior of 
the magnetic Griineisen parameters shown in Fig. [6] 
(g-i), which look qualitatively similar to the ai coeffi- 
cients, but also have some important differences. We 
first discuss the behavior of F.^ for the end compounds. 



The first remarkable feature is the complete absence of 
anomalies in Fq and T^ for C03V2O8 at both Tjjti.i and 
Thtc.i', these F^ coefficients are not constant and both 
strongly increase in magnitude as the temperature is re- 
duced, peaking at the large first-order phase transition 
below which the values decrease again. The behavior of 
Ff, is similar, however a small step is observed in F(, at 
Thti,i- Again, the F^ coefficients of Ni3V208 have oppo- 
site signs as those of C03V2O8, but exhibit qualitatively 
very similar behavior. Also remarkable is the fact that 
both Ff, coefficients have a very similar jump-like fea- 
ture at Thti ■ The absence of an anomaly in F^ at these 
phase transitions implies that the anomalies in ai and Cp 



scale using a single parameter, which however is not con- 
stant but rather increases with decreasing temperature 
(Fig. [6](g-i)). The Griineisen parameters probe the vol- 
ume/pressure dependence of the magnetic energy scales. 
For the simplest magnetic phase transition, one expects 
only the exchange energy to be volume/pressure depen- 
dent, leading to a temperature independent Griineisen 
parameter. Clearly our results show that the situation in 
C03V2O8 and Ni3V208 is much more complicated. Be- 
yond this complexity, the divergent-like behavior of F^ 
is reminiscent of the Griineisen parameters of systems 
close to a quantum phase transition, at which Griineisen 
parameters necessarily divergeji^*^ 



J 



o 
CO 



S ^-1 




(Co, Ni)VO„ 

* ^-x x'3 2 8 

x = 
x = 0.12 
C°3VA ;f = o.20 
x= 28 
X = 0.35 
x=0.50 
x=1.0 




-1 

-2 
20 

10 





10 20 

7(K) 



30 




10 20 

7(K) 



30 



10 20 

7(K) 



FIG. 6: (Color online). Temperature dependence of the relative length {L{T) — L(30K))/L(30K) measured along the a-, b-, 
and c-axes [(a), (b), and (c), respectively], of the thermal expansion coefficient a extracted for the a-, b-, and c-axes [(d), (e), 
and (f), respectively], and of the Griineisen coefficients F extracted for the a-, b-, and c-axes [(g), (h), and (i), respectively], for 
(Goi_:,Ni^)3V208 samples of concentration a; = 0, 0.12, 0.20, 0.28, 0.35, 0.5 and 1. 
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The Griineisen parameters for the intermediate com- types of behaviors are particularly interesting and are 
pounds change, as the ai coefficients, from the Co- to shown in Fig. [T] For x = 0.2, all F^ coefficients appear 
the Ni-like behavior as the Ni content x increases. Two to diverge in magnitude down to Thti,i, where there is a 
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FIG. 7: (Color online). Expanded view on the tem- 
perature dependence of the Griineisen coefficients of (a) 
(Coo.8Nio.2)3V208 and (b) (Coo.65Nio.35)3V208 for the a-, 
h-, and c-axes. 



sharp kink in the curves and below which the magnitude 
decreases again sharply. This behavior is clearly different 
from the pure C03V2O8 case, where the "HTIi" transi- 
tion is hardly visible in F^. However, the overall behav- 
ior looks like what one would obtain if all the transitions 
occurring in Co3V20g were replaced by a single second- 
order transition. For x = 0.35, on the other hand, all 
Ti exhibit a jump- like change at Thti,2, with clear sign 
changes of F.^ occurring for the a- and c-axes. 

All the above features point to the very complex na- 
ture of magnetism in (Coi_a;Nia;)3V208, which results 
from an unusual combination of geometrical frustration, 
a complex set of competing magnetic interactions, possi- 
ble disorder, as well as the doping-induced change of the 
ion moment (from S* = 1 to S* = 3/2) etc. 



4. PHASE DIAGRAM, PRESSURE 
DEPENDENCES, AND DISCUSSION 

Fig. [8] presents the doping - temperature phase dia- 
gram of (Coi_a;Nia;)3V208 extracted from specific heat 
and thermal expansion data (see Section [3] and Appendix 
E]) and from neutron scattering experimentsi^ The dif- 
ferent sets of data agree well and indicate that Thti is 
minimal at a critical concentration Xc — 0.25. The anti- 
ferromagnetic phase is noted HTIi for x < Xc and HTI2 
for X > Xc (with the transition temperatures Thti,i and 



Thti, 2, respectively). At low temperature, a first-order 
phase transition occurs at the critical concentration Xc, 
where the direction of the incommensurate propagation 
vector changes;^ Changes of signs in the thermal expan- 
sion and Griineisen coefficients are observed close to Xc- 
Additionally, an enhancement of the magnetic specific 
heat indicates the presence of strong magnetic fluctu- 
ations at low temperature, which probably result from 
competing magnetic interactions. 

The transition temperatures extracted here for the 
compounds (Coi_a;Ni^)3V208, as well as their order and 
the sign of their uniaxial pressure dependences are sum- 
marized in Table U (see also Appendix |X|. Using the 
Clapeyron relation: 



dT^ AL,y 



AS 



(1) 



where dTx/dpi are the uniaxial pressure dependences of 
a first-order transition temperature T^ , ALi and AS" the 
length and entropy changes, respectively, pi a uniaxial 
pressure applied along i and V the molar volume, we 
calculate the uniaxial pressure dependences of Thtc,i 
and Tfm in C03V2O8 and Tc in Ni3V208 (see Table 
IJ). We find that d\nTc/dp., ~ -bdlnTnTCi/dpi ~ 
—5/2dlnTpM/dpi. The signs of these pressure depen- 
dences indicate that a uniaxial pressure applied along a 
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FIG. 8: (Color online). Magnetic {T,x) phase diagram of 
(Coi_2,Nia;)3V208 extracted from our specific heat (green 
squares) and thermal expansion measurements (black dia- 
monds: a-axis, red circles: fe-axis, and blue triangles: c-axis). 
Full grey stars corresponds to transition temperatures ex- 
tracted from neutron experiments by Qureshi et al.— iS and 
open orange stars to transitions obtained from specific heat 
by Kumarasiri and Lawesi^ Arrows indicate that uniaxial 
pressures applied along a and b are equivalent to Co-doping, 
while a uniaxial pressure applied along c is equivalent to Ni- 
doping. 



or b induces an increase of Thtc,i and TpM and a de- 
crease of Tc, while a uniaxial pressure applied along c 
has the opposite effects. The uniaxial pressure depen- 
dences of a second-order transition temperature can be 
extracted using the Ehrenfest relation: 



dpi ACp 



(2) 



where Aai and ACp are the thermal expansion and spe- 
cific heat variations, respectively, induced at the tran- 



sition. Table HIT] presents the pressure dependences 
extracted for the second-order transitions Thti,i in 
C03V2O8 and Thti,2, Tlti,2, and Tc in Ni3V208. The 
pressure dependences of Thti,i and Tlti,2 have the 
same signs as that of the first-order transitions from the 
same compounds. Thti,2, Tx,2, and Tc can be con- 
sidered anomalous, since their uniaxial pressure depen- 
dences have different signs from that of the first-order 
transition Tc- 



TABLE I: Values of the transitions temperatures, orders of the transitions, and signs of the uniaxial pressure dependences 
along a, b, and c of the transition temperatures. 



Ni-content 

X 



Transition 
temperature 



(K) 



Order 



Sign of 
dn/dpa 



Sign of 
dT./dpb 



Sign of 

on /dp. 




(C03V2O8) 



(NigVaOg 



Thtia 
Thtc,i 
Tlti,i 
Tfm 



ThTI,2 

Tx,2 

TlTI,2 

Tc 
Tc 



11.25 ±0.10 
8.55 ±0.10 
6.2 ±0.1 
6.00 ±0.05 



second 
first 
first 
first 



9.5 ±0.1 
7.20 ±0.05 
6.5 ±0.1 
3.90 ±0.05 
2.45 ±0.10 



second 

second 

second 

first 

second 



± 
± 



± 
+ 

+ 



0.12 


Thtia 
Tl 


8.0 ±0.3 
5.0 ±0.5 


second 
second 


± 
± 


± 

± 


- 


0.20 


Thtia 


7.0 ± 0.5 


second 


± 


± 


- 


0.28 


Thti,2 


6.0 ±0.4 


second 


- 


- 


+ 


0.35 


Thti,2 


6.3 ±0.4 


second 


- 


- 


± 


0.50 


Thti.2 


7.5 ±0.3 


second 


- 


- 


± 



+ 



+ 

+ 
+ 



r 



The magnetic transitions in the mixed compounds are 
not associated with well-defined step-like anomalies in S, 
Li, Cp or ai, and an extraction of their uniaxial pres- 
sure dependences cannot be done as for the parent com- 
pounds. We will limit here the discussion to the signs 
of the uniaxial pressure dependences of Thti (summa- 
rized in Table |T]), which are given by the signs of the 
thermal expansion anomalies at Thti- The trend indi- 
cated by the first-order transitions in the parent com- 
pounds is followed by the doped compounds. Indeed, 
dTHTi,i/dpa,dTHTi,i/dpa > and dTHTi,i/dpc < on 
the Co-side(a; < Xc) while, dTHTi,2/(^Pa,dTHTi,2/dpa < 
and dTHTi.,2/dpc > on the Ni-side(a; > Xc). Since 
Thti.i decreases with x while Thti, 2 increases with x 
(cf. phase diagram of (Coi_2;Nia;)3V208 in Fig. [5]), a 
first approximation is that the effects of a uniaxial pres- 
sure along c are similar to that of Ni-doping, while the 
effects of a uniaxial pressure along a or 5 are equivalent 



to that of Co-doping. This simple picture is compatible 
with the signs of most of the pressure dependences of the 
parent compounds, i.e., for all the first- and second-order 
transitions of C03V2O8, and for the transitions Tlti,2 
and Tc in Ni3V208. Thti,2, T^,2, and Tc in Ni3V208 
are exceptions to the picture proposed here, since their 
pressure dependences are not consistent with it. 

To understand the magnetic properties of 
(Coi_2;Nia;)3V208, a microscopic description of the 
magnetic interactions and of their uniaxial pressure 
dependences is now needed (cf. Yi_j;Laj;Ti03, where 
uniaxial pressure dependences of Tc and T^r were 
related to distortion-induced modifications of the super- 
exchange [29|). However, the high number of parameters 
is such that solving the (Coi_a;Nia;)3V208 problem 
appears to be one of the most challenging puzzles in 
modern magnetism. Indeed, both parents compounds 
are subjects to a cascade of magnetic transitions below 



12 K, for which commensurate and incommensurate 
wavevectors have been reported, indicating the presence 
of a high number of magnetic interactions. Magnetic 
frustration, but also low-dimensional magnetism, are 
key ingredients to understand the magnetic properties 
of this family of compounds. The fact that Ni3V2 08 is 
made of 5* = 1 ions, while C03V2O8 is made of S* = 3/2 
ions, affects surely the magnetic interactions and the 
magnetic single-ion anisotropy. It should be considered 



J 



to understand the properties of the parent compounds, 
as well as that of the mixed compounds. The kagonie- 
staircase lattice of the system is also made of two 
non-equivalent magnetic sites. Finally, a consequence of 
this high number of complex magnetic properties, the 
presence of a multiferroic phase in Ni3V2 08, has still an 
enigmatic origin and will certainly deserve again a lot of 
efforts before being understood. 



TABLE II: Step-like variations of the entropy AS and the relative lengths ALa/La, ALb/Lt, and AL^/Lc at the first-order tran- 
sition temperatures Thtc,i and Tfm in C03V2O8 and Tc in Ni3V208, and uniaxial pressure dependences of these transitions 
extracted using the Clapeyron relation. 



Ni-contont Transition 
X temperature 


AS 
(mJ/molK) 


ALa/La 


ALt/Lt ALc/Lc 


a(ln T^)/dpa 
(kbar-i) 


d\n T^/dpb 
(kbar-i) 


91n T^/dpc 
(kbar-i) 


Thtc,i 
(C03V2O8) Tfm 


0.04 ±0.01 
3.19 ±0.30 


0.004 ± 0.001 
0.54 ±0.07 


0.009 ± 0.002 -0.010 ± 0.002 
0.80 ±0.10 -1.13 ±0.15 


0.11 ±0.05 
0.24 ± 0.05 


0.21 ±0.10 
0.36 ±0.08 


-0.26 ±0.11 
-0.51 ±0.12 


1 (Ni3V20s) Tc 


1.08 ±0.10 


-0.26 ±0.03 


-0.34 ±0.04 0.705 ±0.080 


-0.51 ±0.11 


-0.67 ±0.14 


1.41 ±0.29 



TABLE III: Step-like variations of the specific heat divided by temperature ACpT and of the thermal expansion Aaa, Aat, and 
Aoc at the second-order transition temperatures Thti,i in C03V2O8 and Thti,2, Tlti,2, and Tc' in Ni3V208, and uniaxial 
pressure dependences of these transitions extracted using the Ehrenfest relation. 



Ni-content Transition 


A{Cp/T) Aaa 


Aoi, 


Aac a(ln T^)/dpa 


, 91n T^/dpb 


ain T^/dpc 


X temperature 


(mJ/molK) (10-« R-^) 


(10-« K-i) 


(10-6 K-i) (10-2 kbar-: 


') (10-2 kbar-i) 


(10-2 kbar-i) 


(C03V2O8) Thtia 


480 ± 70 8.24 ± 0.50 


12.8 ±0.8 


-12.9 ±2.0 1.2 ±0.2 


1.6 ±0.3 


-1.6 ±0.4 


ThTI,2 


1005 ±100 1.3 ±0.2 


-9.1 ±0.8 


11.4 ±1.0 0.12 ±0.03 


-0.8 ±0.02 


1.0 ±0.2 


(Ni3V208) TlTI,2 


1080 ±100 -3.35 ±0.50 


-2.25 ±0.25 


6.49 ±0.50 -0.4 ±0.1 


-0.300 ± 0.005 


0.8 ±0.1 


Tc 


881 ±70 10 ± 1 


-2.98 ±0.40 


-9.1 ±1.0 3.8 ±0.7 


-1.1 ±0.2 


-3.5 ±0.4 















CONCLUSION 



In summary, we have examined in detail the 
magnetic phase transitions and phase diagram of 
(Coi_a;Nia;)3V208 usiug heat capacity and high- 
resolution thermal expansion measurements. For the 
parent compounds C03V2O8 and Ni3V208, the same se- 
quences of transitions as that reported in the literature 
are observed. Our heat capacity data show that there are 
no further phase transitions below 2 K. There is a smooth 
evolution of entropy with Ni substitution, going from the 
full S ~ 1 Ni^+ magnetic entropy in Ni3V2 08 to half of 
the S = 3/2 Co^+ magnetic entropy in C03V2O8. The 
thermal expansion data together with the heat capacity 
data allow us to calculate the uniaxial pressure effects 
on the transition temperatures without actually apply- 
ing any real pressure. We find that the uniaxial pressure 
effects for both C03V2O8 and Ni3V208 are extremely 
anisotropic along the three crystallographic directions. 



Intriguingly, the anisotropy of C03V2O8 and Ni3V208 
are just the opposite, suggesting an intimate connection 
between the magnetic ordering in both compounds. This 
is a little surprising, since the incommensurate antiferro- 
magnetic structure in C03V2O8 and Ni3V208 have prop- 
agation vectors along different directions and Co- and Ni- 
ions have different spin states. Detailed calculations of 
the magnetic exchange interactions under uniaxial stress 
are needed to explain this. The thermal expansion data 
suggest that uniaxial pressure should be a good tuning 
parameter, i.e., that applying uniaxial pressure along the 
c-axis of C03V2O8 will lead to a suppression of aU phase 
transitions, whereas c-axis pressure in Ni3V208 is ex- 
pected to increase the transition temperatures. It would 
be particularly interesting to suppress the magnetic or- 
dering in C03V2O8 to zero temperature by c-axis uni- 
axial pressure and then to see if a different magnetic 
order emerges at even higher pressures. Unfortunately, 
one would need about 5 GPa, which is experimentally 



not very feasible. Closely related to the pressure depen- 
dences of the transition temperatures are the magnetic 
Griineisen parameters, which for C03V2O8 and Ni3V208 
exhibit a divergent-like behavior over a limited tempera- 
ture interval reminiscent of a system close to a quantum 
critical instability. The anomalous Griineisen parameters 
are greatly reduced in magnitude for the intermediate 
concentrations, probably due to disorder effects. How- 
ever, a small sign change of the Griineisen parameters is 
observed for intermediate compositions near Xc — 0.25, 
which may be the signal of the remnant doping induced 
transition between the two propagation directions. The 
cascade of magnetic ordering transitions already pointed 
to a quite complex ordering scenario in both C03V2O8 
and Ni3V2 08, and the complex behavior of the magnetic 
Griineisen parameters observed in the present work cor- 
roborates this complexity. 
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Appendix A: Focus on the transitions 

Here, we detail how we have extracted the characteris- 
tic temperatures, as well as the jumps AS", AL^/L^, ACp, 
and Atti, at the different magnetic phase transitions of all 
(Coi_a;Nij;)3V208 samples studied in the present work. 

Figure |9] shows how the different first-order transition 
temperatures have been extracted for each S and Li data 
curve, as well as the jumps AS* and lS.Lil Li at the tran- 
sitions. Figure [10] shows how the different second-order 
transition temperatures of each sample have been ex- 
tracted from each of our Cp and a^ data curve, as well 
as the jumps ACp and Aa^ at the transitions. 

Table HVl presents the values of the different transition 
temperatures extracted here from each data curve (Cp, 
Ua, afc, and Uc). In Table HI an average value (with an 
appropriate error bar) has been retained for each transi- 
tion. 
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FIG. 9: (Color online). Focus on the anomalies at the first-order transitions temperatures in the magnetic entropy and length 
variations. 



10 



o 



-^20 
0.12^ 


HTl,^ 


j 


a) 




l\ V 


/ 


, 


L^r 


X 


K 


/ (Co,V,0,) 


^„.;^^ 




^ 


J x=a2g,^ 

t'> t 

1 V 


H7/,1 ,/ 

y 




7„ l(=°3V.O.) 


\l 




c-axis 


/ 


■^ ■ 




T(K) 



7(K) 



7(K) 



7-(K) 



FIG. 10: (Color online). Focus on the anomalies at the different second-order transitions temperatures in the heat capacity 
and thermal expansion. 

TABLE IV: Transition temperatures extracted here from each data set (specific heat, thermal expansion along a, b and c, and 
nature (first- or second-order) of the transitions. 



Ni-content 


Transition 


Order 


from Cp 


from Qfa 


from Q(, 


from ac 


X 


temperature 




(K) 


(K) 


(K) 


(K) 




Thti,i 


second 


11.26 ±0.02 


11.12 ±0.05 


11.31 ±0.05 


11.38 ±0.05 





Thtc,i 


first 


8.50 ±0.05 


8.50 ±0.05 


8.55 ±0.05 


8.60 ± 0.05 


(C03V2O8) 


Tlti,i 


first 


6.21 ±0.05 


- 


- 


- 




TpM 


first 


6.01 ±0.05 


6.05 ±0.05 


6.01 ±0.05 


6.00 ± 0.05 


0.12 


Thti,i 


second 


8.0 ±0.1 


8.0 ±0.3 


8.05 ± 0.30 


8.0 ±0.3 




Tl 


second 


- 


5.3 ±0.2 


4.45 ±0.50 


5.15 ±0.20 


0.20 


Thti.i 


second 


7.45 ±0.20 


6.9 ±0.3 


6.8 ±0.3 


6.7 ±0.3 


0.28 


ThTI,2 


second 


5.65 ±0.30 


6.2 ±0.3 


6.3 ±0.3 


6.05 ± 0.30 


0.35 


ThTI.2 


second 


6.05 ±0.30 


6.35 ±0.30 


6.28 ±0.30 


6.35 ±0.30 


0.50 


ThTI,2 


second 


7.45 ±0.50 


7.5 ±0.3 


7.5 ±0.3 


7.55 ± 0.30 




ThTI.2 


second 


9.42 ±0.02 


9.54 ±0.05 


9.52 ±0.05 


9.55 ±0.10 




n,2 


second 


- 


7.13 ±0.05 


7.25 ±0.05 


7.16 ±0.05 


1 


TlTI.2 


second 


6.48 ± 0.02 


6.60 ±0.04 


6.50 ±0.05 


6.59 ± 0.05 


(NiaVaOg) 


Tc 


first 


3.93 ± 0.02 


3.93 ±0.05 


3.94 ± 0.05 


3.90 ± 0.05 




Tc 


second 


2.42 ±0.03 


2.53 ±0.05 


2.50 ± 0.05 


2.48 ± 0.05 
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